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Abstract 

A microscopic theory of the plasma resonance in layered metals 
is presented. It is shown that electron-impurity scattering can sup- 
press the plasma resonance in the normal state and sharpen it in the 
superconducting state. 



1 Introduction. 

The plasma resonance can be observed as a steep drop of reflectivity when 
the frequency of an electromagnetic wave goes above the threshold plasma 
frequency Q. The latter is determined in an isotropic plasma or a metal by 
the well-known equation: 

n^ = ^ (1) 

where n is the density of the carriers (electrons), m is their effective mass and 
8 denotes high frequency dielectric constant due to the inner atomic shells. 

If a metal can become superconducting, its London penetration depth A at 
zero temperature is related to the plasma frequency via a simple relationship: 
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where c is the speed of hght. 

Experimental observation of the plasma resonance in conventional met- 
als is difficult because the plasma frequency usually occurs in the far ultra- 
violet range, where inner photoeffect processes are substantial. The situation 
changes drastically for layered metals when the field polarization is perpen- 
dicular to the layers. Indeed, one can naively apply eq. (P with a large 
effective mass M corresponding to the tunneling motion across the layers. 
However, we show that the naive approach fails due to the openness of the 
Fermi-surface. The correct answer is at least m/M times less (m being a 
small mass for in-plane motion). Thus the plasma frequency VLz for c-axis 
polarization occurs in the far infrared or even in the millimeter range. 

Polarized reflection measurements from Y Ba2CuzOj crystals have been 
carried out by Koch et al For the radiation polarized perpendicularly 
to the Cu02 planes a striking difference has been discovered between the 
reflectivity spectra in the normal and in the superconducting state. Typical 
spectra below the superconducting transition (QO-ft') show a sudden sharp 
drop from unity to 0.4 in the vicinity of 12meV (90cm~^) at T = lOK and to 
0.55 in the vicinity of 8meV (GOcm"^^) at T = 65K whereas at temperatures 
above reflectivity just decreases smoothly from 0.8 to 0.60 — 0.65 in the 
whole far infrared region — 15meV. 

A similar behavior of the far infrared c-polarized reflectivity from the 
La2-xSrxCu04^ single crystals has been reported by Tamasaku et al [Q]. 
Plasma edge followed by a dip has been observed by them in the super- 
conducting phase at even lower frequencies (25cm~^ at x = 0.10 and 50cm~^ 
at X = 0.16). This spectral feature being very sharp at low temperature be- 
comes smeared when the temperature is increased and eventually disappears 
in the vicinity of the critical point. 

The physical reason for this phenomenon can be readily understood. 
In the normal state a strong scattering of the electrons by impurities and 
phonons suppresses the plasma oscillations. However, in the superconduct- 
ing state there exists a superfluid of Cooper pairs which are not scattered, 
and thus maintain the coherence of the plasma wave. 

The idea of abnormal transparency of layered superconductors has been 
advocated by T.Mishonov [||. He pursued a phenomenological approach 
based on the two-fluid model of superconductivity (Gorter-Kazimir model). 
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In this model the dielectric function is postulated in the following form: 



The first term in the r.h.s. of this equation describes the contribution of 
the valence electron bands and high frequency optic oscillations. The second 
term is due to the Cooper pairs and the last one is contributed by the Fermi 
excitations. The superconducting and the normal plasma frequencies, Qg and 
fin, are related to the corresponding BCS densities Ug and in the same 
way as in eq. (|I|). The absorption is characterized by some temperature- 
dependent attenuation coefficient 1/r. This approach has been first employed 
for the analysis of the optical properties of the high-Tc superconductors by 
de Marel a/ 0. A recent theoretical work by Tachiki et al goes along 
the same lines. 

The purpose of the present work is to present a microscopic theory of 
the plasma resonance in the normal and the superconducting states of lay- 
ered metals. The anomalously low value of the plasma frequency for the 
c-polarized radiation is a straightforward outcome of this theory. Besides, it 
provides an interesting and important insight into the role played by elastic 
scattering of electrons by impurities. The scattering puts down the refiection 
of electromagnetic waves by destroying the coherence of plasma oscillations 
in the normal state once the scattering rate is comparable with the plasma 
frequency. In the superconducting state, however, the frequent scattering 
prevents the normal carriers to be involved into the collective motion of the 
superconducting Cooper pairs. Therefore, plasma oscillations of the charged 
superfiuid being virtually independent maintain the refiection below plasma 
edge. 

The paper is organized as follows. We introduce the underlying model 
in Section |^. The refiectivity in the normal state is evaluated in Section ^ 
using the kinetic equation. This approach is insufficient to treat the same 
problem in the superconducting state. Not only the quasi-particle occupation 
numbers, but also the Cooper pair wave function varies in the electromagnetic 
wave (see 0,0), so that no local dynamics exists for superconductors. An 
outline of the field-theoretical calculations accounting most naturally for the 
above-mentioned variations and for the random fields of impurities is drawn 
in Section H. The results of these calculations are represented and discussed 
in Section ^. We return to the discussion of advantages and limitations of 




(3) 
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our model and summarize its physical implications in the last Section ^. The 
analytic continuation of the electromagnetic response function is carried out 
in [Appendix A| . The electromagnetic response for the frequencies above the 



threshold 2 A is studied in [Appendix B 



Preliminary version of this work has been published in the Proceedings 
of SPIE conference. 



2 The model. 

Most layered conductors have a rather complicated crystal structure includ- 
ing, in particular, several conducting layers {CUO2) per elementary cell. 
These layers are usually coupled much more strongly than those belonging 
to different cells. In this article a simplified model with a single layer per 
elementary cell is employed. To make the model even simpler, the electron 
dispersion relation is taken to be totally isotropic in the a6-plane. Although 
the spectrum of the real compound possesses only orthorhombic symmetry, 
the anisotropy in the a6-plane is not very pronounced. In other cases the 
symmetry may be tetragonal. That gives us hope that the isotropic model 
provides a reasonable approximation. 

In our model the electrons are assumed to be scattered predominantly 
by the impurities. The collisions are assumed to be elastic, and since the 
overlapping of the electronic wave functions even between the neighboring 
planes is small, only in-plane scattering is taken into account. In principle, 
two competing types of processes, that of in-plane and inter-plane scattering 
with different collision times r and r', respectively, should be considered 
independently. Note that the inter-plane collision time r' must be compared 
to the effective time of inter-plane quantum tunnel hopping suppressed by 
the in-plane scattering [Q. This time is equal to (7^t)~^ where 7 is an inter- 
plane hopping amplitude. We assume the strong inequality r' 3> (7^t)~^ 
to be valid, allowing us to disregard the inter-plane scattering. Our theory 
relies on the assumption that the metallic properties of the CuO planes are 
not spoiled by scattering, i.e. e^r ^ 1 {ep stands for the Fermi energy). On 
the other hand, the dimensionless parameter 7r measuring the relative value 
of the inter-plane tunneling with respect to the in-plane scattering may vary 
in a broad range. The role of inelastic processes and other more realistic 
modifications to the theory will be considered elsewhere. 
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To summarize, our model is that of non-interacting massive electrons 
propagating freely in the ab-plane and hopping to the neighboring planes. 
Their motion in the ideal crystal is governed by the following dispersion law 

(^{p,Pz) = p'^/2m + 7 (1 - cos (pzd/h)) (4) 

where p = {px,Py) is an in- plane component of momentum and Pz is a quasi- 
momentum component parallel to c-axis. The effective mass and the inter- 
layer spacing are denoted by m and d respectively. Finally 7 is the hopping 
amplitude introduced earlier. 

The following obvious relationship between the hopping amplitude 7 and 
the effective mass in c-direction M 

7 = h^Md'' (5) 

will be employed later. 

As far as the superconducting properties are concerned only the simplest 
ground state characterized by an isotropic and homogeneous distribution 
of the order parameter is considered. The BCS approximation is adopted. 
Although the potential drawbacks of the a weak coupling theory for the 
quantative description of high-Tc superconductors are clearly understood, 
we believe that the basic features of the plasma resonance phenomenon are 
revealed in the framework of this approach. 



3 Reflectivity in the normal metal. 

We now evaluate the reflectivity for the incident wave propagating along 
X-axis, the surface of the layered metal coinciding with the yz-plane, and 
electric field being normal to the layers, i.e. directed along 2;-axis (c-axis). 
The current Jz{uJ, k^) induced by the electric field Ez{u}, k^) in the skin layer 
is also directed along the 2;-axis: 

Jz{u;, k:c) = (Tzz{'^)Ez{tO, k^,) (6) 

where uj and k^ are the frequency and the wave number of the incident wave, 
respectively. The conductivity (Tzz depends in general both on wave number 
and frequency 

( i.\ o- 2 [ dn^^^ v^jpjd^p , . 

(Tzz{uJ, kx) = - 2ze / — — — 7 
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where n^^\e) is the Fermi occupation number. Neglecting space dispersion 
{kvp compared to a; + i/r ) one obtains a kind of Drude-Lorenz formula 

^Uu;) = . (8) 

ivr/i [uo + i/t] 

Let us note that r in the above expression is a true collision time in contrast 
with the formula for the homogeneous metals where it is substituted by a 
transport timeQ 

Henceforth we will omit tensor and vector indices of the conductivity 
tensor azz, electric field E^, magnetic field By and wave vector and kx used 
in the present work and denote them simply ex, E, B and fc, respectively. 

The dielectric permeability of the metal is related to the conductivity as 

£{u) = £ + ^^a{u) =£(l- ^} \ (9) 

where E denotes the high-frequency dielectric constant (a contribution due to 
the valence electrons and high-frequency optical modes), plasma frequency 
is defined as 

2e^m , , 

The dielectric permeability, in addition to being a function of cj, depends on 
two additional parameters /3 = vpjc and fi^r. Both of them are intrinsic 
characteristics of the metal. Although /? is always small, fi^r may vary in 
different compounds. We assume that it cannot be too large, obeying the 
following inequality 

/?fi,r < 1 (11) 

If (pJ] ) is satisfied, the spatial dispersion of the conductivity a may be ne- 
glected, as we have done already. Indeed, one can verify that a strong in- 
equality kvp ^1 cj + z/r I is vahd where 



c 



k = -JS{u) (12) 



provided (|TTD holds. 



■^The reason for this difference is the canceUation of the "arrival" term in the colhsion 
integral due to the independence of the scattering amplitudes on pz . 
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The reflectivity is related to the dielectric constant by the equation 

2 



R{iu) 



S{iu) - 1 



S(uj) + 1 



(13) 



Several typical spectra of the reflectivity are shown in Fig.l. The steep drop 
of reflectivity appearing in the vicinity of the plasma frequency Qz when the 
latter is large with respect to the scattering rate 1/r is usually followed by a 
dip 0. The bigger is the high-frequency dielectric constant S, the deeper is the 
minimum. When the impurity elastic scattering rate 1/r is significant, with 
respect to Qz , the corresponding reflectivity spectrum becomes featureless. 

Let us emphasize two peculiarities of the interaction between the electro- 
magnetic field and layered metals. First, the plasma frequency is surprisingly 
small. In fact, eq. (|10|) for may be rewritten as 

where a is an in-plane lattice spacing. The above expression for the plasma 
frequency contains two small factors, m/M and {h/ppd)'^, reducing the con- 
ventional values by two or three orders. Therefore, it is unlikely that fi^r ^ 1 
in experiments on the copper oxides, at least until single crystals of very high 
quality become available. It looks even less likely that the criterion ([TT|) can 
be violated. Hence the Drude expression (|^) may serve as a good approxi- 
mation at any frequency. Recall that the additional small factors appear as 
a consequence of the openness of the Fermi- surface. 

The second peculiarity concerns the absence of the so-called anomalous 
skin-effect for z-polarized radiation. In this anomalous regime, a leading role 
is played by a strong spatial dispersion. The latter might become essential 
once the mean free path of the electrons exceeds the skin-layer depth. How- 
ever, in this case the inequality pflr <C 1 would be violated, whereas it is 
definitely satisfied in the experiments with high-Tc superconductors. 

To conclude this section, we argue that all our results are robust with 
respect to the replacement of the quadratic dispersion e(p) = /2m by some 
arbitrary anisotropic dispersion e(p), provided the Fermi line (in the ah- 
plane) is sufficiently smooth. Indeed, the essential relationship which allows 

''We do not consider 8 < 1. 
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to neglect spatial dispersion is the inequality kvpT 1, which is equivalent to 
the condition (pTl). The integration along the Fermi-line in (|^) gives basically 
the same formula as for the quadratic dispersion. Only the plasma frequency 
is redefined: 

where N{e) is the density of states per unit volume. 

4 Field theory approach. 

For the reason given in the Introduction, we apply field theory methods to 
the calculation of the dielectric permeability of layered superconductors. If 
we set the gap in the spectrum, A, to zero, we get the result for the normal 
state. This justification of the kinetic equation by more rigorous field theory 
methods is useful since the validity of the kinetic equation is not obvious 
when 7r < 1. 

The electromagnetic response tensor Qais relates the average current J to 
the vector-potential A 

Jo,{uj,r) = - J Qo^piu, r, P)Ap{uj, ?)£r' (15) 

The quantity Qapi^-^,^,^') is usually called the Pippard kernel. According 
to the Kubo formula, the response function (Pippard kernel) Qai3i^^,r,r') is 
equal to the Fourier transform of the retarded Green function of two current 
densities: 

^ iT 

Q^p{u;,f,r') = — dte"^\[Ut,r),MO,r')]). (16) 

C JO 

where T is the temperature of the system and the angular brackets denote 
the averaging over the Gibbs ensemble. Following the conventional routine, 
we first compute the auxiliary Matsubara response 

iT r^/T 

Ql'p{uj^,r,r') = - / dTe'--^{UT,r)MO,r')) (17) 
c Jo 

where Um = 2nmT. 

The retarded response (|16]) can be found via analytic continuation of the 
Matsubara response ( [T7| ) from the discrete points u = iuj^ in the upper 
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half-plane of the complex variable cu, using the well-known relation 

Q^f^iiu^, f, P) = Qt'pium, , r, rO (18) 

If elastic impurity scattering only is taken into account, the Matsubara Green 
function of currents in (|17D may be expressed in terms of two-point electronic 
Green functions. We prefer to write down the corresponding relation in mo- 
mentum space rather than in the coordinate space; this will prove convenient 
for subsequent averaging over impurity configurations: 

tI ji^^ji^^''''^'^'''^'^^^^^^ ^^^^ 

where r]^ = {2u + 1)itT, and p± =p±k/2,q± = q±k'/2,P± = P±{k-k')/2, 
and the group velocity Va{p) = d^/dpa- Finally, the vertex F, defined as 

F(u;„,pV,p_,g+,g_) = T^G'^^(u;„ + r]^,p+, q+)G'^^ {r]^,p_, q.) + 

T J2F''{uj.^ + r],,p^,q^)F''{r],,p^,q_) (20) 

depends on the four momenta p±,q±- Let us emphasize that the electronic 
Green functions are the exact Green functions for the non-interacting elec- 
trons moving in the potential of some arbitrary but fixed impurity configu- 
ration, and obeying the well-known Dyson integral equation ||^. Averaging 
over the electronic Fermi distribution has already been performed in ([T9|) , 
but averaging over the random ensemble of impurities is yet to be done using 
the standard diagrammatic technique p|. After this second averaging, trans- 
lational invariance is restored, and the averaged Green functions depend on 
twice as few variables: 

G''{v.,P,q} = S{p-q)G''{v,aP)) 

F^'{Vu,P,q) = S{p + q)F^'{v,m) (21) 

These reduced Green functions were calculated by Klemm et al [|10[, by 
a natural modification of the method of averaging originally proposed by 
Abrikosov and Gor'kov for isotropic superconductors: 
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.Mt^ ^^ _ ir/g(r/)/e(r/) +^ 
e(r/)2 + ^2 



Here 

e(r/) = V'^^ + e~(r7) = e(r/) + l/(2r) (23) 

The averaging of the vertex (pOf ) normally involves a solution of the Bethe- 
Salpeter equation or, more accurately, a system of equations for different ver- 
tices. For the dynamic processes the number of components increases since 
the vertex depends on two independent frequencies and can be retarded or 
advanced with respect to each of them 0, §]. Fortunately, these difficulties 
can be avoided for the component Q^z of the response tensor. The calcula- 
tions in this case are simplified due to the independence of the amplitudes of 
the electron scattering by the impurities on the 2;-component of the electronic 
momentum. In the framework of this approximation for the Bethe-Salpiter 
diagrams, integration over the z-component of the momentum in each of the 
vertices can be separated and involves the following integrals: 

IT /d rjin 

-Tv/d {ZTI ) 

where ^± = ^ ± {k, v)/2. These vanish | because Vz is the total derivative of 
a periodic function of Pz- According to the above arguments the component 
of the response tensor Qzz is given by: 

QfzM = -—[ 7^^,^(p.)r^(^„,e+,e-) (24) 



(2vr) 
where 

+ F^\u^ + r^,,i+)F^\r^,,i^)\ (25) 

To obtain the Pippard kernel, one must perform the analytic continuation 
of the quantity F and then integrate it over S,- The details of this procedure 
may be found in the Appendix A. 

^This nullification corresponds to the vanishing of the "arrival" term in the collision 
integral of kinetic equation. 
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5 Plasma resonance in a superconductor. 

For the specific problem under consideration we neglect spatial dispersion, i.e. 
the dependence of Qzz on the wave-vector k. For the dielectric permeability 
we use eqs. with the following modification: 

47rc 



£{uj) = £- 



-Qz 



(26) 



In the superconducting state the penetration depth does not increase to infin- 
ity at zero frequency due to the Meissner effect, but it still has its maximum 
at zero frequency. We have argued in Section 3 that it is much bigger than the 
mean free path at high frequencies. This condition is satisfied without doubt 
at low frequencies. The following calculations are similar to the treatment of 



the high-frequency response of superconductors by Bardeen and Mattis [12 



and by Abrikosov, Gor'kov and Khalatnikov . We modified the method 
of the latter authors to account for the layered structure and the impurity 
scattering. 

After the analytical continuation in the [Appendix A| (see eq. ( [A.?] )) of the 
vertex T^{uj) in the eq. (^4]), the integration over ^ is performed explicitly 
in the eq. ( |A.13| ). We express the Pippard kernel Qzz 

e^7^m d 



Qz 



(27) 



via the dimensionless function K^{uj, T). The formulae for its real and imag- 
inary part in the frequency range uj < 2A derived straightforwardly from the 
eq. ( [A.13| ) are presented below: 

TleK«(.,T) = r tanh(../2rM„ [. , + A'^ , r 

Ja-^ el + (e + 1 tY 



£6 1 



+ : 



+ l/r 



I fOO 

XmK^{uj,T) = — dr] 
2r J A 



- 1 




e + - 



tanh ^ + tanh 

2T 2T 



1 V 
tanh — 

2T 



dr] (28) 




(tanh^ -tanh^ ) (29) 
V 2T 2TJ ^ ' 
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where the notations has been changed as follows 



r]+ = r] + uj e = - A2 e+ = \/r]l - e = ^A^ - rj^ (30) 

Although the normal limit of the Pippard kernel at finite frequency cannot 
be extracted from the above relations, we have verified (see discussion on the 
page 22 following the eq. ( |B.2| )) that 



A=o {uj + i/r) 



(31) 



confirming our earlier calculation based on kinetic equation. 

Let us analyze the frequency behavior of the kernel in several regimes cor- 
responding to different values of temperature and different levels of scattering 
rate. The main difference between the superconducting and the normal states 
is that ^22(0;) vanishes for u; = in the normal state, but it takes a positive 
finite value in the superconducting state in conformity with the London rela- 
tion between the current and the vector-potential. The London penetration 
depth A2 is connected to the value of the kernel Qzz{0): 

X;' = A7rQU0)/c (32) 

Two first terms in the eq. (|28D for TZeK^{uj) contribute to the static limit 
of the Pippard kernel: 

->'^---4-7:7m^^i?T7i;^ ,33, 

This London factor vanishes as A — > 0. 

The Pippard kernel may be evaluated exactly in the clean limit : 

lim K^(uj,T) = 1 (34) 

r >oo ^ ' 

Notice that the limiting value doesn't depend on temperature. This result 
is naturally associated with the suppression of dissipation. The Cherenkov 
absorption of photons is forbidden since the Fermi velocity is much less than 
the speed of light. On the other hand breaking of Cooper pairs is kinemat- 
ically permitted but its rate is proportional to the Bogolyubov coherence 

factor [T^\ (u^_^_^v^ — v^_^i:Upj where k is the momentum of the photon. For 
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the normal skin effect (A; = 0), which we only consider, this factor is equal 
to zero. 

Comparing the clean limit to the static one we come to conclusion that 
in a weak scattering regime (Ar ^ 1) and for the temperature relatively 
close to the transition point the real part of the Pippard kernel grows rapidly 
from its static value (p3D to unity. This onset becomes less pronounced at 
lower temperatures (Fig 2A). The increase of scattering rate is accompanied 
by reduction of the mean level of TZeK^uj) and smothering of its variation in 
the low frequency range (Fig 2B). 

At higher frequencies there is a plateau K^{uj,T) ^ 1 ending in the 
vicinity of the threshold u = 2 A. This is the point of a mild logarithmic 
singularity 

^ 7rr(ix;-2A), /2A-u;\ , A , . , 

TZeK [u, T) = log ( — — — I tanh — + regular terms (35) 

where the derivative with respect to u becomes logarithmically large. This 
singularity is closely related to the Cooper pair breaking above the threshold 



[see discussion on the page |2^) 



Thus, by changing the temperature one can observe a large variation of 
the plasma frequency when 2A approaches the bare plasma frequency ujp 
but still is larger than the latter. The attenuation can remain reasonably 
small to provide a sharp plasma resonance since Q, r and Tc are independent 
parameters. 

The observable plasma frequency at a given temperature is defined as a 
root of equation 

ne£{uOp) = (36) 

In a sufficiently clean superconductor (A r ^ 1) and for the moderate values 
of the bare plasma frequency Q (|14D the observable plasma frequency ujp 
almost coincides with the bare one. However, ii Qt < 1 the observable 
plasma frequency Up diminishes reaching its static limit at -C 1/r. In 
more dirty samples the deviation of the observable plasma frequency from the 
bare one is very significant in the whole range < Q < 2 A. The diminution 
of temperature causes the augmentation of the observable plasma frequency. 
Nevertheless, in a very dirty limit (A r -C 1) the observable plasma frequency 
is by a factor y/ At smaller than the bare one even at zero temperature. 
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The behavior of the real part of the Pippard kernel is illustrated in Fig.2A 
(r Tc = 0.1) and Fig.2B (rTc = 3.4). Two plots are presented on both fig- 
ures: dashed for T/T^ = 0.3 and dashed-dotted T /Tc = 0.74 respectively. 
Looking at the dirtier case one can see that the observable plasma frequency 
is less than its static limit c/A at the lower temperature and exceeds this 
limit at the higher temperature. It would be very instructive to study exper- 
imentally the temperature variation of the plasma frequency in more detail 
and to compare it to the temperature dependence of the London penetration 
depth in c-direction A^. The latter may be obtained independently from the 
measurements of the radio-frequency impedance. Our theory predicts ex- 
tremely weak temperature dependence of the observable plasma frequency in 
clean samples whereas this dependence may be fairly remarkable in samples 
with a moderate concentration of impurities. 

The imaginary part of the Pippard kernel is determined by the dissipa- 
tion. As it has been already mentioned, two types of photon absorption 
are substantial: the Cherenkov absorption of a photon by a quasiparticle 
(Landau damping) and the absorption accompanied by the dissociation of a 
Cooper pair. The former process is forbidden kinematically in the absence 
of scattering, the latter is suppressed for the photon momentum k = 0. At a 
finite scattering rate the momentum is not conserved any more, and both pro- 
cesses become possible. The pair dissociation requires the photon frequency 
u! > 2 A. This range of frequencies is considered in [Appendix B| . 

The probability of the Cherenkov absorption is proportional to the Fermi 
occupation numbers of quasiparticles. These vanish exponentially at T — > 0. 
The asymptotic form of Xm for A ^ T and u ^ T can be readily found 



from eq. (pOf) : 



T ^ -A/T • r ^ W f ^\ 2Ar^2 (2 + ^/A) 

where Kq is a MacDonald function of the zeroth order. Examining it one can 
find that the dissipation depends substantially on the scattering rate and in 
a somewhat surprising way. Namely, it becomes especially small provided 
Ar ^ 1. Even as the temperature approaches to the transition point and 
the number of the excitations increases the dissipation may continue to be 
small if the scattering rate is sufficiently high. The physical reason for this 
phenomenon, as already mentioned in the Introduction, is the decoupling of 
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a strongly dampen excitation motion from the collective plasma oscillations 
of the Cooper pairs. However, the asymptotic formula (0) does not apply 
in this case. 

Several plots of imaginary part of the Pippard kernel versus frequency are 
shown in Fig. 3. A minimum of XmK^{uj) occurs at some frequency which 
depends on Tc, T and r. Several characteristic frequency scales for XmK^{uj) 
can be easily seen in eq. (^). These scales are T, A, and (Ar^)~^. When 
T is the smallest of the scales the minimum of ImK^iuj) is located at ~ T 
and the minimal value of ImK^{uj) may reach the order of unity for the 
moderate temperatures in the clean case. We also see that ImK^{uj) is 
presumably a linear function of a; at -C T and behaves as tu"^/^ at S> T, 
but still uj < A. One can readily find numerous types of limiting behavior of 
XmK^{uj) depending on the relationship between above-mentioned scales. 

Reflectivity spectra are shown in Fig.4. The plots illustrate the behavior 
of the spectra with respect to the temperature variation and the dependence 
on a dimensionless parameter^ rTc, intrinsic for a specific sample. The values 
of two other sample-specific parameters: high frequency dielectric constant 
and bare plasma frequency are chosen to be £^ = 25.0 and fi^/Tc = 1.2 
for each of the four plots. Each plot corresponding to some scattering rate 
(rTc = 0.1; 0.9; 1.7; 3.4) is combined from three subplots: a solid curve for 
the normal metal and two for the superconducting state at temperatures 
T /Tc = 0.3 and T/Tc = 0.74. At the largest value of the scattering rate 
(Fig.4A) plasma edges followed by the "beaks" are clearly visible on both 
subplots below while the normal subplot is completely featureless. Subse- 
quent plots show the evolution of the plasma feature for the higher tempera- 
ture. At the lowest scattering rate tTc = 3.4 (Fig.4D) the normal spectrum 
is still smeared whereas the superconducting spectra for both temperatures: 
T = 0.74 Tc and T = 0.3 Tc are already close to the limiting coUisionless 
spectrum. 

6 Conclusions. 

We have shown that the plasma resonance can be suppressed in the normal 
state of a layered metal due to a strong electron-impurity scattering, but it 
should be restored in the superconducting state since the superconducting 

^We choose the transition temperature as a natural energy scale. 
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part of the electron liquid is not subject to scattering. Moreover, below 
scattering sharpens the plasma edge because it reduces the attenuation by 
removing Fermi-excitations from the coherent wave motion. 

This feature is common for both the BCS and the two-fluid (Gorter- 
Kazimir) models. Indeed, the imaginary part of the dielectric function (|^) 
is small for the resonance frequency ujp = Qg provided a strong inequality 
fl'^T -C fig is satisfied. This condition is obviously violated in the vicinity of 
the superconducting transition. 

Let us further discuss the relationship between the present theory and the 
two-fluid model. Namely, we are going to address the question, to what extent 
the phenomenological expression @ may approximate the analytic structure 
of the real dielectric function which we believe is correctly described by eqs. 



(p6|-p9|). This issue is important since the Kramers- Kronig transformation 
used routinely in the analysis of the experimental reflectivity data is very 
sensitive to the singularities of the analytic function involved. 

Let us first compare the behavior of the real part of the dielectric func- 
tion in a weak scattering regime: At ^1. In this case the BCS behavior 
described in the previous section is reasonably well reproduced in the frame- 
works of the two-fluid model. In particular, the static limit of uj'^TZeS{uj) is 
proportional to the density of superconducting Cooper pairs. The growth of 
the response real part due to the contribution of the normal carriers saturates 
at the same level as in the BCS model when a; S> 1/r. However, even in the 
clean limit the shapes and (which may be more important) the characteris- 
tic scales of the low frequency onset are different. Indeed, the intermediate 
frequency scale 1/A^r arising naturally in the BCS theory cannot be found 
in the Gorter-Kazimir model. Finally, it is impossible to incorporate the 
threshold phenomena in the vicinity of = 2A, like a logarithmic singu- 
larity of TZeK^{uj,T) described by the eq. (|35|), into the two-fluid plasma 
electro dynamics . 

The augmentation of scattering worsens the agreement between the two 
models. In a strong coupling regime: Ar <C 1 the two-fluid model predicts 
an anomalously long plateau at the level of the static limit because the nor- 
mal electrons are overdampen and do contribute only to the imaginary part 
of the Pippard kernel. The BCS strong coupling behavior of the TZeK^{uj) 
is different (see Fig. 2B). In fact, the two- fluid model encounters some prin- 
cipal difficulties at small enough Ar. It is not clear how the normal and 
superfluid densities can be deflned consistently when the non- conservation of 
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the momentum becomes significant. 

In the work the clean limit of the two-fluid model has been considered. 
This approach looks dubious because the fact that the resonance is suppressed 
in the normal state combined with the experimental value of uOp ~ 2/3 A 
imply that Ar < 1. 

As far as ImK^{(jj) is concerned one can find less common features than 
the discrepancies between the two-fluid and the BCS models. It is not sur- 
prising, since the imaginary part of the Pippard kernel as a function of uj 
has only one characteristic scale, r~^, whereas it depends on four different 
scales: r~^. A, T, (Ar^)~^ in the present theory. In particular, the position of 
the minimum of the function XmK^{uo) in the two- fluid theory it is always 
^min = In the BCS calculations it depends strongly on the temper- 

ature. At low temperature ujmin ~ T and is always smaller than r^^. A 
striking difference between the BCS and the two-fluid dissipation in a strong 
scattering regime (rTc = 0.1) is demonstrated in Fig. 5A. However, even in 
the clean superconductor the dissipation calculated in the frameworks of the 
BCS model is much weaker than its two- fluid counterpart (Fig 5B). Of course 
the light absorption through the Cooper pair breaking remains beyond the 
scope of the two-fluid model. 

Summarizing, the two-fluid model gives a reasonable, though in some 
range of variables rather crude approximation for TZeS{uj) and fails to de- 
scribe correctly Xm£{uo). 

Although many realistic features are lacking in the underlying model we 
believe that our main findings will survive in its more refined and sophis- 
ticated versions. It is likely that the in-plane anisotropy of the electronic 
dispersion would change the macroscopic constants, but not the shape of the 
plasma resonance. It is not so clear, however, what would be the impact of 
the anisotropy of the superconducting gap. One may conjecture e.g. that if 
the d-wave pairing mechanism of superconductivity prevails the attenuation 
at low temperatures would not be exponentially small. 

In a more ambitious treatment, one may have to reexamine the scat- 



tering model. According to the experimental data []TS| the resistance in 
the c-direction displays a semi-conducting rather than a metallic behavior. 
Therefore, the inter-plane motion might be better described in terms of hop- 
ping between the defects than between the planes. We argue, however, that 
at high frequencies long-distance hopping is ineffective, whereas short-range 
hopping is regulated by inter-plane tunneling. Nevertheless a detailed treat- 
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ment of these problems is highly desirable. 

Although the plasma resonance has been observed only in the supercon- 
ducting state so far, future developments in a technology of single crystal 
preparation could drastically raise the chances of its observation in the nor- 
mal state. Such single crystals could be utilized as 100% polarizers in the far 
infra-red region. 

It would be very interesting to search for the plasma resonance in the 
more pronounced layered superconductors, such as Tl- and Bi-based oxides, 
where the edge could occur in the microwave region. Especially promising 
would be epitaxially grown multilayers, in which the transition temperature 
and the energy gap are determined by the superconducting layers whereas 
the plasma frequency is determined by the square of the inter-plane tunneling 
amplitude. 
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discussions. Our special thank is due to 1. Bozovic and W. Saslow for their 
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Appendix A Analytical continuation. 

According to Section | (eq. (0)), to calculate the Pippard kernel one must 
perform an analytic continuation of the quantity F defined by eq. (^) which 
is reproduced here for visibility omitting the factor —2ie^/c. 

+ T^F^^K + r/.,e+)F^^(r/.,e-) (A.l) 

from a discrete set of points lOn = i2mTT with n > into the whole upper 
half-plane of the complex variable uj. We have restored the spatial dispersion 
omitted in the text, since the analytic continuation can be performed in a 
more general context. 

Once analytic continuation is performed the integration over ^ is in turn. 
Let us note, however, that the integral over ^ taken in the infinite limits and 
the sum over i/ together are divergent. Instead of introducing an explicit 
cutoff in course of calculations one can deal with convergent expressions 
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employing a proper regularization procedure. To this end we subtract from 
T{uj) the same quantity a.t u = and A = 0. The corresponding Green 
functions will be labeled by an additional subscript 0. 

We start from the analytical continuation of the Green functions calcu- 
lated by Klemm et al (|2^ , |23|) from a discrete set of arguments rj^ = i7r(2z/+l). 
First, we rewrite these functions in the following form: 

(6(r^) + l/2rf + e 



(e(r/) + l/2rr + e2 

Expanding the above expressions over the poles of the formal variable e{r)) 
and analytically continuing the latter as a function of the variable t] one 
obtains the retarded Green functions in the following form: 



2 \eR{r]) + i/2r - ^ e^(r/) + i/2T + 
= 27m {eR{v)+i/2T-^ + e^(r^)+V2r + e) ^^'^^ 

where 

eR{r]) = v'(r7 + iO)2-A2 (A.4) 

The branch of the square root is chosen to give positive values of eji{ri) for 
7] > A. The infinitesimal imaginary number is added to its argument in order 
to prevent the branching points and the poles of the above functions from 
entering the upper half plain of the complex variable rj. The advanced Green 
functions can be obtained from the retarded one by a formal substitution 
r — >• — r and en^r]) — > 6^(77) where 

eAiv) = ^{V - - A2 (A.5) 

According to the above consideration in Section 4 the component of the 
response tensor Qzz is given by the following expression 

Q^K,fc) = (^^.'(p.)rK,e+,e-) (A.6) 
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where f 2 = {2ti'-) / d) s\.n{2TT'-)Pz/ d). 

Keeping the external frequency in the subset of discrete imaginary num- 
bers: uj = iuJn = i2TmT we can transform the sum in the r.h.s. of eq. ( |A.1| ) 
into an integral: 



2T 



drj + 



^ J^^ G^(r/ + ioon, e-) tanh d^ + 

_L ^^(^ + i^^^ e-) tanh d ^ 



(A.7) 



where the contours Fi, in the complex plane ri go around the intervals 

of the imaginary axis from = iO to = icxo, from r/ = —iiOn to r/ = iO 
and from 77 = —zoo to = — ici;„ respectively, all counterclockwise. These 
contour integrals can be readily transformed into the integrals along the real 
axis, using the periodicity of tanh(?7/2T) with the period uJn- The last step 
is to put + iO instead of luJn in all integrals to get the "retarded" analogue 
of the vertex ( [A.l|) written in the following integral form 



drj 
~oo 47ri 



tanh — - 
2T 



\2T 



tanh 

VZIJ 



\2T 



F«^(u;,r^,e+,e-) 



(A.8) 



where double-retarded, double-advanced and mixed vertices are defined by 
the following equations: 



(A.9) 

(A.IO) 
(A.ll) 



Here and in what it follows ri± = r] ± uj. The expressions (|A.9 - A.lTl) are 
obviously analytic in the upper half-plane and since all the integrals are now 
convergent, one can safely perform the integration of it over ^ first. This 
integration becomes especially easy provided the representation ( |A.3| ) and 
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its analogue for advanced Green functions are used. The contribution due to 
the subtracted terms vanishes since 



1 



T] + i/2T - i 
1 



(A. 12) 



The result of integration of vertex r^(ti;, (see eq. ( |A.8|) ) over ^ depends 

on the scalar product kv. The remaining integration over (p and Pz becomes 
trivial, once the spatial dispersion is neglected, leading to the following final 
expression: 

r/+r/ + \ tanh(?/+/2T) - tanh(r//2T) 



dr] 

-oo 4 



— 1 



+ tanh {^^ 



ejj(r7)eK(r?+; 



eA{v)^AiV- 



(A.13) 



The real and imaginary parts of K^{uj) for u < 2 A are given by eqs. 
of the text. See Appendix B for the expressions of the same values and their 
analysis in the range u > 2 A. 



Appendix B Above the Threshold. 



Here we present formal expressions for the real and imaginary parts of the 
normalized Pippard kernel in the range u > 2A and a brief analysis of their 
behavior at zero temperature. The required expressions can be extracted 
straightforwardly from the general equation ( |A.13| ): 



neK{u 




rj+T] + A^ 



tanh ^ + 
2T 



tanh — + tanh — 
2T 2T 



2T 



tanh 



2T 



dr] 
T 



+ 



+ 
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/A y 



2 e + e_ 



+ 



l/r2 



tanh — - + 
2T 



fl - !^!L±^) ftanhA _ tanh^^l (B.l 
l/r2 \^ ee_ y V 2r 2T J \ ^ 



1mK{uj 



tanh(r//2r) 
e-e_)2 + l/r2 
A2 




(^tanh ^ - tanh ^ 



+ 



tanh ■ 



2T 



(B.2) 



The notations are the same as in the text. 

Let us first take a look at the normal limit A — > 0. In this case the first 
integral in the l.h.s. of the eq. (p.lj ) just vanishes whereas the coherence 
factors (1 ± {rjrj± + /\^)/ee±) are reduced either to or to 2 due to the 
following identities: e+ = r^^, e = and e„ = — r/-. The Lorentzian factors 
associated with the non- vanishing coherence factors tend to ujjiuj^ + I/t^) 
or (l/r)/(u;2 + l/r^) in formulas for the real and the imaginary parts of the 
kernel respectively. Taking into account the above remarks one can evaluate 
the remaining integrals straightforwardly and finally arrive to the Drude- 
Lorentz formula (0). 

We specially consider the case of zero temperature. Then equations ( [B.l| ) 
and ( [B.2|) look much simpler: 



neK{u) 




rj+rj + A^ 



1 + 



e_|_e 
?7_?7 



dr] + 



A2 



rj^r] + A^ 



drj 



(B.3) 
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1 PLU-A 



dri (B.4) 



The imaginary part of ( p.4|) corresponds to the absorption of a photon ac- 
companied by the creation of two quasiparticles (dissociation of a Cooper 
pair). The probabihty of this process, suppressed in the absence of scatter- 
ers, is proportional to the scattering rate 1/r. In addition to the coherence 
factor it contains the Lorentzian distribution function for the total momen- 
tum of a pair of quasiparticles located on the opposite points of a diameter. 
When r — i> oo, this distribution turns into (5-function. 

The integral in the above formula (|B.4|) for the imaginary part may be 
expressed in terms of the standard elliptic integrals as follows: 



(CU + 2A) (l-4A2r2 + u;2^2) \^ (l + u;2^2 
(u; + 2A)r 4 A r (A + u; + u;^ r2^ 



„ Eimiuj)) + ^ '-.Kimiuj)) (B.5) 

1 + ^ ^ " (c^ + 2A) (l+a;2r2)' ^ ^ ^ ' 

Here K (m), E (m) and 11 (n, m) are complete elliptic integrals of the first, the 
second and the third kind respectively. The moduli of the integrals depend 
on the frequency and the scattering rate. 

/c^-2A\2 ^ ^ r2 (cu - 2A)' (1 + cu2 r2) ^ 

- (^) ' ^) - - (l + (.2-4A2).2) (B-6) 

The function ( p.5|) is always negative. Just above the threshold its absolute 
value increases linearly with the intercept being inversely proportional to the 
scattering rate: 

lmK{uj) = -7rr(cu-2A)/2 + 0(cj/A-2)^ (B.7) 

The non-linear phase of the descending slope may be approximated by means 
of the following interpolating function: 

ImK{uj) ^ I , - 1 I (B.8) 

4Ar L/l + 4(.;-2A) Ar2 ' 
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This interpolation gives a remarkably good estimate of the actual 2mK{uj) 
almost up to the point Umin ~ 2A + 0.9/r where the imaginary part reaches 
its minimum Q A minimal value of XmK{uj) significantly depends on the 
scattering rate. If the latter is small compared to the threshold frequency 
{At ^ 1) the dissipation is also small: XmKiojmin) ~ — 7r/4Ar. However, 
in the regime of a strong scattering the dissipation growth is saturated at 
the level: ^ 

\im^ImK{uj^in) = -- (B.9) 

Examining eq. ( |B.8|) one can easily see the appearance of the intermediate 
frequency scale scale 1/Ar^. 

In the high frequency range the imaginary part of the Pippard kernel 
decays monotonically with the asymptotics: 

ImK{uj) = -1/uj T + 0(1 /uj^) (B.IO) 

in conformity with the Drude-Lorentz formula for the normal metal (|3T|) . The 
existence of the gap becomes irrelevant for the plasma oscillations driven 
by the electromagnetic wave provided the frequency of the latter is large 
compared to the threshold value. 

The discontinuity of the derivative with respect to frequency of the imag- 
inary part of the Pippard kernel at = 2 A implies that the total Pippard 
kernel, considered as a function of the complex variable u, has a logarithmic- 
type singularity at this point. The behavior of this function in the vicinity 
of its singularity is governed by the eq. (P^D. The real part of the kernel 
decreases rapidly when approaching the threshold with a negative logarith- 
mically divergent derivative. 

Switching the temperature on causes the augmentation of the dissipation. 
On the contrary, lZeK{uj) at finite temperature is diminished. Still near the 
threshold the asymptotic behavior ( |B.8| ) is valid if one considers this contri- 
bution as complementary to the Cherenkov absorption part, and multiplies 
it by the factor tanh(A/2T). 
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Figure Captions. 

Figure 1 

Reflectivity versus frequency (measured in relative units: uj/VLz ) in the nor- 
mal state. 

Figure lA: The high-frequency dielectric constant 8 — 25. 

Figure IB: The high-frequency dielectric constant £ — A. 

The different curves on both plots correspond to the different values of scat- 
tering rate Vt^T = 0.01; 0.1; 0.5; 1.0; 5.0. 

Figure 2 

Real part of the normalized Pippard kernel versus frequency (measured in 
relative units: uj/T^ ). The bare plasma frequency fi^ = 1.2 Tg. The high- 
frequency dielectric constant £^ = 25 . Three spectra are plotted for each of 
the following scattering rates: 

Figure 2 A: tT^ = 0.1. 

Figure 2B: tT^ = 3.4 

The dashed plot corresponds to T/Tg = 0.30, A/Tc = 1.73 and the dashed- 
dotted one corresponds to T/T^ = 0.74 , A/% = 1.38 respectively. The 
values of A(T/Tc) have been calculated in the frameworks of the BCS theory. 

Figure 3 

Imaginary part of the normalized Pippard kernel. The bare plasma frequency 
Qz = i-^Tc- The high-frequency dielectric constant £^ = 25 , The dashed 
plot corresponds to T/T^ = 0.30, A/Tc = 1.73 and the dashed-dotted one 
corresponds to T/T^ = 0.74 , A/T^ = 1.38. 

Figure 3 A: tT^ = 0.1. 



Figure 3B: tT^ = 3.4 
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Figure 4 

Reflectivity versus frequency (measured in relative units: uo/Tc ). The bare 
plasma frequency = 1.2 T^. High-frequency dielectric constant £ = 2b . 
Three spectra arc plotted for each of the following scattering rates: 

Figure 4A: tT^ = 0.1. 

Figure 4B: tT^ = 0.9 

Figure 4C: rT^ = 1.7 

Figure 4D: tT^ = 3.4 

The solid curve is a plot for the normal state. The dashed plot corresponds 
to T/Tc = 0.30, A/Tc = 1.73 and the dashed-dotted one corresponds to 
T/Te = 0.74 , A/Tc = 1.38 respectively. 

Figure 5 Comparison of the imaginary parts of the normalized Pippard 

kernel for the BCS and the two-fluid models. The solid curves correspond to 

the BCS model. The dashed plots correspond to the two-fluid model. 

The bare plasma frequency = 1.2 Tc- High-frequency dielectric constant 

£: = 25 , t/t; = 0.30, a/t; = 1.73 . 

Figure 5 A: tT^ = 0.1. 
Figure 5B: tT^ = 3.4 
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